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Getting the position of the instrument in starting traverse and 
staking out surveying can be very helpful for the surveyors. 
The most common method is the placement of the instrument 
on the known point, then those surveys are possible to be 
accomplished. This research is aimed to develop a new 
method and procedure to get x, y, and z values of the 
unknown position of the instrument based on two known 
points and the law of cosines. The method of this research is 
the implementation of the law of cosines and Euclidean 
Norm in solving the problem of getting the coordinate of 
instrument position. The innovation of this procedure has not 
been used yet in survey practice and has not been 
accommodated in electronic distance measuring (EDM) 
based survey instruments such as Total Station. The 
experiment of measurement to test the procedure is 
conducted virtually using the total station of SimusurveyX 
1.0.7. The total measurement of ten random triangles is 60 
times, where each triangle is measured 6 times. The result of 
measurement is close to the ground truth, and it can be 
repeatable. The implication of this research is enabling the 
surveyors to shortcut traverse measurement by locating Total 
Station in the first unknown point of the traverse. 
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1. Introduction 


The cosine law can be explored deeply to get the invention of new methods and breakthroughs in 
determining and predicting the location and position of certain objects. The law of cosine has 
been proved in determining and predicting position, it has been successfully implemented in the 
Coordinate Rotation Digital Computer (CORDIC) Processor for location awareness system and 
robot indoor location technology [1,2]. There are a lot of demands for breakthroughs and new 
practical methods in determining the position of objects for various kinds of applications in any 
subject. Such breakthrough is needed urgently in mapping and surveying subjects to facilitate 
engineers and surveyors in the field to increase their productivity. The breakthrough can be an 
invention of work instruction that can solve the field problems or procedures which is can be 
accommodated in the firmware of mapping and surveying instruments or application in 
ubiquitous devices possible to control the instruments [3]. Some breakthroughs in mapping and 
surveying such as mobile and computer applications have been developed by software engineers 
to enhance the productivity of surveyors in the field [4]. The method to stake out or transfer 
location from the map to the field using are available on the total station software features, and 
there is an available android app to facilitate such measurement using theodolite or low-cost 
instruments [5]. Even recently, the method of learning for surveying subjects using a virtual 
learning environment to enhance student competencies has been developed, and the development 
computer-based tutorials are being developed [6]. 


The famous method in starting terrestrial mapping using a Total Station is setting up the 
instrument in the known position points or known as Benchmark (BM) and setting Back Sight 
(BS) to the known point, or in other words to determine the location of the target point the 
instrument must be set on the known point and refers (BS) to the other known points [7]. Based 
on the operation procedure of Total Station in the mode both of traverse measurement and stake 
out measurement the surveyor must set up the total station instrument on the known position 
point and direct the telescope to the known BS point, and the next surveyor can determine Fore 
Sight (FS) point or locate Stake Out (SO) point. The other method determines roughly the 
position or coordinates of the point from a couple of known positions by the trilateration method 
[8,9]. A similar method also has been simulated in MATLAB to get the coordinate of an 
unknown point using the three-point resection problem [10]. 


Most wide range products of total stations facilitate users to determine the coordinate of target 
points based on known backsight reference and the known location of the instrument. The 
research question in this study is how do we know the current random location of the instrument 
if there are two available known point positions, and how such kind of method can shortcut the 
traverse’s measurement in certain field conditions. The aim and objective of this research are to 
develop a new procedure for determining the position of the Total Station in x and y coordinates 
based on two known position reference points using the law of cosines and how to apply its 
method in traverse measurement. The implementation of the method of instrument location 
determination in practical traverse surveying needs accuracy and precision. The accuracy can be 
defined as the closeness between the agreement or accepted reference value or the truth to a large 
number of the test results or observations, meanwhile, precision is the closeness of the test 
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results or consistencies of the observation groups and evaluated based on size differences [11— 
14]. 


This study proposes a new method of determining position x, y, and z coordinates of random 
unknown position points based on the distances and angles of two known position points. Such a 
method has not been accommodated in the method of the stake-out or traverse survey. The 
Impact of this newly proposed method can solve the field problem when the common procedure 
of BS cannot be implemented due to the field obstacle. In the future, this new procedure should 
be accommodated by the latest survey instruments in the built-in application to enhance their 
performance and will facilitate their users in the field. This new method will enrich the 
conventional method of the stake-out and traverse survey using a common instruments such as 
Total Station and Theodolite. 


3. Methods 


The SimuSurveyX version 1.0.7 is hired in this research activity, in the simulation of 
measurement experiments, the software is developed by the V-Lab (Visualization Laboratory) 
Department of Civil Engineering National Taiwan University and can be downloaded freely at 
https://simusurveyx.caece.net/software/. SimusurveyX provides a realistic effect and learning 
environment [15]. In running the simulator and the app needs a PC or Laptop with Windows XP; 
Windows Vista or Windows 7/10 operating system and Minimum system requirement : 1.6 GHz 
CPU, 384 Mb RAM, 1280 x 720 display. The method of this research consisted of installing 
the simulator, designing measurement scenarios, calculating the position of the instrument, 
validating results, determining accuracy and precision, and discussing and drawing conclusions. 
This method of development is based on the law of cosines which can be explained as follow: 


Based on Figure 1, if A, B, and C are the Angles and a, b, and c are the length of sides, therefore 
it is governed by the following formula [1] which the formula has been successfully 
implemented in Coordinate Rotation Digital Computer in Location Awareness System. 


c? = a* + b* — 2ab Cos C (1) 
b* = a* + c* — 2ac Cos B (2) 
a* = b* +.c* — 2be Cos A (3) 


The length and angle of the triangle can be solved using Formula 1, 2, and 3 [16]. 


Based on the law of cosines theorem, the Total Station or Instrument Position can be determined 
precisely through the following procedure of measurement: 


1) Setting up the total station in P3 at any location that it has no barrier Line of Sight to P1 
as known position and P2 as known position (See Figure 4) 

2) Aligning the telescope to P1 as Backsight (BS) set horizontal angle 0° 0’ 0” measure 
horizontal distance between P3 and P1 note as b (Refers to Figure 4). 
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3) Rotate telescope to P2, note horizontal angle as C, and measure the distance between P3 
and P2 as a, 

4) Calculate the length of c or horizontal distance between Pl and P2 using triangular 
principle or the Pythagorean theorem [2] or Euclidean Norm [17] that can be written as 
formula (4) 


C= V(X, — %1)* + (V2 — 1)? (4) 


5) Based on the law of cosine [1] formula (2), Calculate the angle of B 


a?+c?—b? 
B=arc COR ae (5) 


6) Calculate the azimuth of P1 to P2 


a = arctan ar (6) 
if X2_ x, > Oand yz — y; > 0 then azimuth (d1_,p2) = @ (7) 
if X2_ x, > Oand yz — y, < 0 then azimuth (0p1_,p2)= 180-4 (8) 
if X2_ x, < Oand yz — y, < 0 then azimuth (dp1_,p2)= 180 + a (9) 
if X2_ x, < Oand yz — y; > 0 then azimuth (01 _,p2) = 360 — a (10) 


7) Calculate azimuth P2 to P3 
Ap23p3 = 4pisp2 — B + 180 (11) 


8) Determine the Coordinate of P3 (x3, y3, z3) using this formula 


X3 = X_ + aSiN Ap2-,p3 (12) 
¥3 = V2 +t ACOS Ap2,p3 (13) 
Z3 = Z_ — (atan(90 — V) + Hi — Hr) (14), 


i.where V is vertical Angle P1 to P3, Hi is the height of the instrument, and Hr is the height of 
the reflector or mirror 

9) Validate the Coordinate, the result can be validated using the position of Instrument in 
SimuSurveyX by clicking the escape key and then selecting the result in the screen menu. 

10) Determine the error of x and y by comparing mathematical calculations and the 
instrument’s location determined by the simulator. 

11) Repeat the experiment of measurement 

12) Determine the accuracy and precision of the coordinate of P3 determination using this 
method. 


1 
Faccuracy vt ie (Ri — Rerue) (15) 
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Precision = R + Ina eg; — BR)? (16) 


13) Assessing accuracy is also can be done using Root Mean Square Error (RMSE) [18] and 
Mean Average Error using Formula 17 and 18 


1 2 

RMSE = |- Pa (%pi — You) (17) 
1 

MAE == Yina|Ypi — You (18) 

P1 

(x1, Yi) 
P2 
ar? ¢ : (X2, Y2) 


Fig. 1. Two Known Points and The Unknown Point of Total Station Position. 


Such calculation simply can be accomplished using a spreadsheet formula, or with little effort, it 
can be solved by developing certain software by transforming the procedures into algorithms, 
and algorithms into code. In this research as a tool to get the result from random experiments, the 
Microsoft Excel formula is developed to solve the problem and compared the result using a 
simple computer application. 


4. Results 


From the instrument setting in the simulator, the position of the mirrors and instrument is 
depicted in figure 2. The x, y, and z of point A (P1) or mirror 1, point B (P2) or mirror 2 and C 
(P3) or station 


of the instrument is ground truth. 


The default error of the virtual instrument axis is set to 0.000 m, as depicted in figure 3, where 
the x-axis is 0.000 m, the y-axis is 0.000m and the z-axis is 0.000 m. 
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ERROR SETTING 


| X-Axis: 0.00 
Y-Axis: 0.00 

| Z-Axis: 0.00 
RESET 3-AxISs ERROR 
DISTANCE ERROR: ON 


Back 


lei 


Fig. 4. (a) Mirror 1, and (b) mirror 2 reading. 


In this measurement experiment base on the simulator, the location of rulers and instruments are 
depicted in figure 4. Where Ruler 1 at point A or P1 and Ruler 2 at point B or P2 


xX, = 6.880 m X) =17.640 m 
yi=9.640 m y2=17.440 m 
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Z1=30.010 m Z2=30.129 m 
meanwhile, figure 4 (b) can be known as following 
C=51° 41’ 17” = 51.688° 
The distance of P3 (instrument) to P2 (mirror 2) 
Instrument height (Hi) = 1.267m 
mirror 2 height (Hr) = 1.168 m 
a= 10.545 m 
the distance of P3 (instrument) to P| (mirror 1) 
mirror | height (Hr) = 1.164m 
b= 16.936 m 


length of c or the distance of P1 to P2 can be calculated using formula 4 


c = (17.640 — 6.880)? + (17.440 — 9.640)? = 13.290 m 


and using the law of cosine formula 5, the angle of B can be calculated 


10.5452 + 13.290 * — 16.9367 


R= 
Te COS > x 10.545 x 13.290 


= 89.798° 
The azimuth of P1 to P2 (ap; yp2) can be calculated using formula 6 up to formula 10 as follow 


17.640 — 6.880 
17.440 — 9.640 


= 54.601° 


a =arctan 


Where the condition full fills x,;_ x2 > 0 and y, — yz > 0, then the Azimuth P1 to P2 is 
Opi _,p2 = N 54.601° E 
The next step is, the calculation of the azimuth P2 to P3 (ap2_,p3) using formula 11 
Op1_,p2 = 54.601 — 89.778 + 180 = N144.263° E 
The Coordinate of P3 (x3, y3, z3) using this formula 12 and 13 
x3 = 17.640 + 10.545 sin 144.263 = 23.799m 
y3 = 17.440 + 10.545 cos 144.263 = 8.881 m 
Z3 = Zz — (10.545 tan(90 — 92.499) + 1.267 — 1.168) = 30.490 m 
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The above calculation can be solved manually using a calculator, spreadsheet, or computer 
application. The proposed procedure is transformed into Programming Code and deployed as a 
simple computer application. The enhancement of functionality and performance of the 
application is an improvement and is needed to further develop different platforms such as 
android and IOS [19]. Both calculations using a spreadsheet and developed application have the 
same result, the development of the application can be adopted by the Total Station manufacturer 
so that the user or surveyor can run this procedure directly from the instrument. The simple 


developed computer application to calculate the new procedure is depicted in figure 4. 


Table 1 


ng Intrument Position 


Finding Unknown Posiition of Instrumen (Total Station) 


Based on Two Known Point 


Known Position 


P1 Coordinate: 


Height of rmument (Hi) 


P2 Coordinate 


@ 1706 


gm irae 


Height of Retector (Hr) 


1.108 


Horizontal Distance 
aP3toPZ | 10045 


qpramPs (Was 


Cuca | 


i= 


Fig. 4. Computer Program to solve the procedure that is developed using Java. 


Comparison of x, y, z Measurement, and Ground Truth. 


# PA Ground Truth P2¢round Truth P3 Measurement P3 Ground Truth Error (P3y-P3¢r) 
xl yl zl x2 y2 Z2 x3 y3 z3 x3 y3 z3 x3 y3 z3 
6.880 | 9.640 |30.010 |17.640 | 17.440 |30.129 |23.799 | 8.881 |30.490 | 23.80} 8.880 | 30.487 |-0.001| 0.001 | 0.003 
6.880 | 9.640 |30.010 |17.640 | 17.440 |30.129 |23.802 | 8.884 |30.491 | 23.80] 8.880 | 30.487 | 0.002] 0.004 | 0.004 

1 6.880 | 9.640 |30.010 |17.640 | 17.440 |30.129 |23.799 | 8.882 |30.490 | 23.80} 8.880 | 30.487 |-0.001| 0.002 | 0.003 
6.880 | 9.640 |30.010 |17.640 | 17.440 |30.129 |23.802 | 8.879 |30.491 | 23.80] 8.880 | 30.487 | 0.002 | -0.001 | 0.004 
6.880 | 9.640 |30.010 |17.640 | 17.440 |30.129 |23.801 | 8.881 |30.490 | 23.80} 8.880 | 30.487 | 0.001} 0.001 | 0.003 
6.880 | 9.640 |30.010 |17.640 | 17.440 |30.129 |23.802 | 8.879 |30.490 | 23.80} 8.880 | 30.487 | 0.002 | -0.001 | 0.003 
7.257 | 4.860 |30.020 |28.333 | 20.804 |30.239 |22.252 | 3.911 |30.229 | 22.25} 3.912 | 30.225 |-0.002| -0.001 | 0.004 
7.257 | 4.860 [30.020 [28.333 | 20.804 [30.239 [22.255 | 3.915 [30.228 | 22.25] 3.912 | 30.225 | 0.001] 0.003 | 0.003 

2 7.257 | 4.860 [30.020 [28.333 | 20.804 [30.239 [22.255 | 3.912 [30.229 | 22.25] 3.912 | 30.225 | 0.001] 0.000 | 0.004 
7.257 | 4.860 |30.020 [28.333 | 20.804 |30.239 |22.255 | 3.914 |30.229 | 22.25] 3.912 | 30.225 | 0.001] 0.002 | 0.004 
7.257 | 4.860 [30.020 [28.333 | 20.804 [30.239 [22.256 | 3.912 [30.229 | 22.25] 3.912 | 30.225 | 0.002] 0.000 | 0.004 
7.257 | 4.860 [30.020 [28.333 | 20.804 [30.239 [22.253 | 3.912 [30.230 | 22.25] 3.912 | 30.225 +-0.001| 0.000 | 0.005 
9.088 | 12.118 [30.020 [25.742 | 20.336 [30.122 [18.655 | 2.707 [30.840 | 18.65] 2.706 | 30.836 | 0.003] 0.001 | 0.004 

3 9.088 | 12.118 [30.020 [25.742 | 20.336 [30.122 [18.650 | 2.705 [30.841 | 18.65] 2.706 | 30.836 |-0.002| -0.001 | 0.005 
9.088 | 12.118 |30.020 |25.742 | 20.336 |30.122 {18.654 | 2.706 |30.842 | 18.65] 2.706 | 30.836 | 0.002] 0.000 | 0.006 
9.088 | 12.118 [30.020 |25.742 | 20.336 [30.122 [18.652 | 2.707 [30.840 | 18.65| 2.706 | 30.836 | 0.000] 0.001 | 0.004 
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30.050 


0.003 


-0.001 


0.002 


11.35 


20.054 


30.050 


24.620 


11.377 


30.626 


9.802 


2.086 


30.052 


9.802 


2.088 


30.050 


0.000 


-0.002 


0.002 


11.35 


20.054 


30.050 


24.620 


11.377 


30.626 


9.801 


2.087 


30.054 


9.802 


2.088 


30.050 


-0.001 


-0.001 


0.004 


11.35 


20.054 


30.050 


24.620 


11.377 


30.626 


9.804 


2.087 


30.053 


9.802 


2.088 


30.050 


0.002 


-0.001 


0.003 


8.091 


17.286 


30.010 


29.721 


3.851 


30.289 


14.114 


4.946 


30.245 


14.11 


4.943 


30.244 


-0.002 


0.003 


0.001 


8.091 


17.286 


30.010 


29.721 


3.851 


30.289 


14.116 


4.945 


30.244 


14.11 


4.943 


30.244 


0.000 


0.002 


0.000 


1} 8.091 


17.286 


30.010 


29.721 


3.851 


30.289 


14.114 


4.943 


30.245 


14.11 


4.943 


30.244 


-0.002 


0.000 


0.001 


8.091 


17.286 


30.010 


29.721 


3.851 


30.289 


14.116 


4.941 


30.245 


14.11 


4.943 


30.244 


0.000 


-0.002 


0.001 


8.091 


17.286 


30.010 


29.721 


3.851 


30.289 


14.116 


4.942 


30.245 


14.11 


4.943 


30.244 


0.000 


-0.001 


0.001 


8.091 


17.286 


30.010 


29.721 


3.851 


30.289 


14.117 


4.944 


30.245 


14.11 


4.943 


30.244 


0.001 


0.001 


0.001 
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The list of java language code for button action event functions is as follows: 


private void btCalculActionPerformed(java.awt.event.ActionEvent evt) { 


// TODO add your handling code here: 


double xl = Double.parseDouble (txl.getText()); 
double yl = Double.parseDouble(tyl.getText()); 
double zl = Double.parseDouble(tzl.getText()); 
double x2 = Double.parseDouble (tx2.getText()); 
double y2 = Double.parseDouble (ty2.getText()); 
double z2 = Double.parseDouble (tz2.getText()); 
double hr = Double.parseDouble (thr.getText()); 
double hi = Double.parseDouble(thi.getText()); 


double a = Double.parseDouble(ta.getText()); 


double 


oO 
ll 


Double.parseDouble(tb.getText ()); 


double vD = Double.parseDouble (tvD.getText()); 


double vM = Double.parseDouble(tvM.getText()); 


double vS = Double.parseDouble(tvS.getText()); 


double V = vD + (vM/60)+ (vS/3600); 


double azimuth12=0; 


double azimuth23=0; 

double c = pow((pow(x2 - x1,2))+(pow(y2 -yl1,2)),0.5); 
double deg = (pow(a,2)+pow(c,2)-pow(b,2))/(2*a*c); 
double B = Math.toDegrees (Math.acos (deg) ); 


double azl = (x2 - x1)/(y2 - yl); 


double az=Math.toDegrees (Math.atan(azl1l)); 


//double azimuth12 = az; 


double x3 = 0; 
double y3 = 0; 
double 23; 


if ((x2 - x1)>0 && (y2 - yl)>0) { 


azimuth12 


aZ; 


azimuth23 azimuth12 - B +180; 


x3 = x2 + a* (Math.sin(Math.toRadians (azimuth23))); 
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y3 = y2 + a*(Math.cos (Math.toRadians (azimuth23) )); 


else if ((x2 - x1)>0 && (y2 - yl1)<0){ 


azimuth12 180 - az; 


azimuth23 


azimuth12 - B +180; 
x3 = x2 + a* (Math.sin(Math.toRadians (azimuth23) )); 
y3 = y2 + a*(Math.cos (Math.toRadians (azimuth23) )); 
jelse if ((x2 - x1)<O && (y2 - yl)<0) { 


azimuth12 = 180 + az; 


azimuth23 azimuth12 - B + 180; 
x3 = x2 + a* (Math.sin(Math.toRadians (azimuth23))); 
y3 = y2 + a*(Math.cos (Math.toRadians (azimuth23) )); 


selse if ((x2 - x1)<0O && (y2 - yl)>0) { 


azimuth12 360 - az; 


azimuth23 


azimuth12 - B + 180; 
x3 = x2 + a*(Math.sin(Math.toRadians (azimuth23))); 


y3 = y2 + a*(Math.cos (Math.toRadians (azimuth23) )); 


z3 = z2 - (a*Math.tan(Math.toRadians(90 - V))+(hi - hr)); 
DecimalFormat df = new DecimalFormat ("###.###") ; 
tc.setText (String.valueOf (df.format(c))); 

tx3.setText (String.valueOf (df.format (x3))); 


ty3.setText (String.valueOf (df.format (y3))); 


tz3.setText (String.valueOf (df.format(z3))); 
} 


And then check the Error of measurement and calculation by comparing x3, y3, and z3 of 
measurement to the ground truth of P3 to get the error value. It is known station P3 has x3 = 
23.800 m, y3 = 8.880 m, and z3 = 30.487 m. The error of x; measurement is 0.001 m, the error of 
y3 is 0.001 m, and the error of z3 is 0.003 m. In this research 10 triangles with two different 
known locations of points and random total station positions were set, each triangle is measured 
6 times to determine accuracy and precision. Table 1 is the results of the measurement, then 
calculated the difference in x or dx is the result of x3 -x2, dy, and dz, mean, standard deviation, 
and relative standard deviation as listed in table 2. 


The Accuracy and precision in the experiment of measurement 1: 
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Average Accuracy, = 0.001 m 
Average Accuracyy = 0.001 m 


Average Accuracy, = 0.003 m 


SD, = 0.001 m therefore, precisions of dx measurement is 6.161 + 0.001 m 
SDy = 0.001 m therefore, precisions of dy measurement is 6.161 + 0.001 m 


SD, = 0.001 m therefore, precisions of dz measurement is 6.161 + 0.001 m 


Tabel 2 
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Value of Average Errors and % Error of x y z compared to the true differential distance value of (P3 — P2). 


Measurement (i) Measurement (meter) Average P3y— P3¢r 
AXp2->p3 dy p2->p3 dz p2->p3 x y Z 

Average 6.161 -8.559 0.361 0.001 0.001 0.003 
1 SD 0.001 0.002 0.000 

Average -6.079 -16.891 -0.010 0.000 0.001 0.004 
2 SD 0.001 0.001 0.001 

Average -7.089 -17.629 0.719 0.001 0.001 0.005 
3 SD 0.002 0.001 0.001 

Average -18.013 -6.434 -0.948 -0.001 0.001 0.002 
4 SD 0.001 0.002 0.001 

Average 12.898 14.301 0.507 0.000 -0.002 0.005 
5 SD 0.001 0.001 0.001 

Average 19.004 5.610 -0.286 -0.001 0.000 0.004 
6 SD 0.002 0.001 0.001 

Average 9.102 -11.999 0.293 0.000 -0.002 0.003 
7 SD 0.002 0.002 0.001 

Average -8.594 -15.110 0.718 0.000 0.000 0.005 
8 SD 0.001 0.001 0.000 

Average -14.818 -9.289 -0.573 0.000 0.000 0.003 
9 SD 0.002 0.002 0.001 

Average -15.605 1.092 -0.044 0.000 0.000 0.001 
10 SD 0.001 0.002 0.000 

Grand Average 0.000 0.001 0.003 


The accuracy and precision of the measurement are graphically depicted in figure 5 where the 
figure 5 (a) x, y scatter plot shows dots, close to each other in the center of x = 0.000 and y 
=0.000, while 2 other scatter plots (b) and (c) where vertical axis represent z, dots close each 
other but the center in z = +0.004, (d) in three dimensional scatter view showing several dots 
occupy on the above of spherical graphs (radius=0.005m, and center x=0.000m, y=0.000m, 
z=0.000m) means inaccurate, Meanwhile (e) all dots close each other and covered by sphere 
graph which depicts the measurement precision and the center of the assemblage of dots are 
(0.000, 0.000, 0.004) means z measurement less accurate than x and y. 


The accuracy of measurement using the newly proposed method has proof that horizontally has 
high accuracy + | mm and the result of measurement repeatable with a maximum standard 
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deviation of 2 mm. The less accurate issue in z value is a common problem for the instrument 
such as Total Station, the vertical error depends on the vertical angle, the greater the vertical 
angle then the result a higher error in z measurement [20], the higher accuracy is needed the 
height difference measurement can be determined using auto-level [21]. 


In the common procedure of the SO and traverse measurement where the instrument of Total 
Station must occupy the known position point and reading mirror or prism to the other known 
position point as reference or BS, then rotate the telescope to the unknown position point as 
target point that will be determined its position. In BS measurement the Total Station needs the 
free line of sight (LOS) so that the EDM transmits the beam to the prism, in which the prism 
returns the beam to receiving optic and will be converted to an electrical signal [7]. When field 
obstacles such as a tree, or building are found, that causes a barrier between the first benchmark 
point and the second benchmark known point, so the instrument and prism have no LOS. In such 
a field problem, the new proposed procedure can be applied as a solution. The implementation of 
the new procedure in surveying is depicted in Fig. 6. 


Scatter Plot of x, y Error in Meter Scatter Plot of x, 2 Error in Meter Scatter Plot of y, 2 Error in Meter 
0015 } 0015 0015 


901} 00k0}—_ ooie-+— 


Fig. 5. The 3D scatter Plot of x, y, and z error in meter (a) shows horizontally x and y precise and 
accurate, meanwhile (b, c, d) vertically z less accurate and but (e) relatively precise. The sphere has a 
radius of 0.005 meters. 


The accuracy and precision also can be assessed from the scatter plot measured location of P3 
(x3, y3, z3) versus the Ground Truth of P3 (x3, y3, z3) as illustrated in Figure 6. It shows the 
scatter dots relatively fit the line of accuracy and the position of dots in each group of 
observations shows the result is replicable which is prove the precision of the new method in 
general. The value of RMSE and MAE from measured and Ground Truth is calculated based on 
the formulas 17 and 18 using the R Studio function as listed in Table 3. This shows the most 
accurate measurement is on the y-axis while the value of the position on the z-axis is the most 
inaccurate, although in general, the measured value has high accuracy which is proven by the 
small value of RMSE and MAE. 
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x3 Measured vs. X3 Ground Truth 


y3 Measured vs. y3 Ground Truth 
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Fig. 6. Plots of Measured value versus Ground Truth. 
Table 3 
Value of RMSE and MAE of Measured x3, y3, and z3 result of the function (rmse) and function (mae) in 
R Studio. 
Position Coordinate RMSE MAE 
X3 0.001678293 0.001383333 
y3 0.001839384 0.001383333 
Z3 0.003796929 0.00345 
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Traverse without obstacle in between Bm1 & Bm2 


; i 
A i Pe ae eee a P3 
Bm1 el - 
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Bm2 \ ya 
= — 
P4 
P5 


Traverse without obstacle in between Bm1 & Bm2 
(applying new proposed method) 


i we 9 P3 
obstacle Bm2 a it 
A Known Point | 
@® Traverse Point P5 P4 


Fig. 7. A common method to measure traverse measurement [7] in the upper, and the lower is the 
proposed new method that can be used when there is an obstacle between two known points. 


5. Discussion 


The result of a virtual measurement can be calculated in finding the predicted location, 
meanwhile, the simulator has accurate information of ground truth where the instrument is set in 
the field. Such a feature of SimusurveyX has better information than the real grid field, the 
simulator will detect every movement of the instrument with millimeters accuracy, and such 
detection will be reflected in the northing, easting, and elevation coordinates of the instrument 
station. Therefore, those features can be utilized to prove the proposed new method in the self- 
defining location where it will be a little bit more difficult when must be proven using the real 
instrument in the real field. The weakness of this newly proposed method is inaccuracy in 
vertical or elevation measurement, it could be caused by the inaccuracy in determining the height 
of the reflector, as the reflector height is not automatically informed by the simulator, but it must 
be determined using virtual ruler, the virtual ruler is place side by side to the reflector, and then 
reflector height is defined based on the scale of the virtual ruler. The finding of the new method 
can be translated into the algorithm and the algorithm can be translated into the certain 
application. In testing the proposed new method, simple software based on java has been 
developed. The software development in the self-positioning system using the newly proposed 
method has a great opportunity to enrich the standard built-in application software in Total 
Station because there is no availability of built-in software features of this newly proposed 
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method in the recent product of Total Station. There is a great opportunity for users of Total 
Station with the Android Operating System, the user is easier to develop their app to enhance the 
functionality of the instrument. Users can create android apps in a situation where they find 
problems or obstacles in starting the traverse procedure and data processing. In replacing 
standard spreadsheets and desktop applications those have more hardship to use directly in the 
field. 


This proposed new procedure is also can be implemented when surveyors use another instrument 
such as theodolite in the field. The strengths of the new proposed procedure gives higher 
accuracy and precision horizontally it also will reduce time in setting the instrument in the field, 
as the instrument does not need to be set right on the BMs or pegs. And the weakness of the new 
procedure is the vertical position has lesser accuracy rather than the vertical position. So this new 
method will be matched when the surveyor just needs x and y coordinates as depicted in Fig. 5. 
According to the strength and weaknesses of the proposed new procedure, it is suitable for 
planimetric and cadastral measurements. When users need higher accuracy in height differences 
or elevation, they can combine the measurement using Total Station for x and y coordinates and 
auto-level for z coordinates. In the future, this procedure can be used to develop a positioning 
system instrument based on distance and angle from two known positions. Its possible 
application is not merely just in the surveying engineering field, but it has a chance to be applied 
in the navigation system and other problems in a self-positioning system. 


6. Conclusion 


Both traverse and stake-out surveys sometimes find new obstacles in conducting setting the 
instrument. The common obstacle there is a barrier between two known points so the instrument 
does not have a line of sight to the backsight point. Such conditions need a new procedure. The 
new proposed procedure in determining position based on two known position points and the law 
of cosine has proven to determine the coordinate of instrument position, the procedure is tested 
using virtual measurement data and it has high accuracy in x and y between 0 to 2 mm, 
meanwhile vertically z has lower accuracy between | to 5 mm. This new procedure can be used 
to overcome the problem in stake out and polygon measurement when the surveyor finds an 
obstacle between two known points and the old procedure cannot be applied in the field. The 
new procedure starts by setting Total Station or Theodolite in a random unknown position as 
long as the instrument can read prisms or rulers in two known points. Then the location of the 
instrument could be determined accurately and precisely. 
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